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The methods of the probability theory have been used in order to build up a new model of 
hysteresis which is different from the well-known Preisach model. It is assumed that the system 
consists of large number of abstract particles in which the variation of an external control parameter 
(e. g., the magnetic field) may result in transitions between two states and The state 

of a particle is characterized by the value +1 or —1 of a random variable (e. g., the magnetization 
direction parallel or antiparallel to the magnetic field). The transitions are governed by two further 
random variables corresponding to the S^~' =>■ <S' + ' and the S' + ' =>• S^' transitions (e. g., 
"up switching" and "down switching magnetic field"). The method presented here makes possible 
to calculate the probability distribution and consequently the expectation value of the number of 
particles in the (or iS 1 - - - 1 ) state for both increasing and decreasing parameter values, i. e., the 
hysteresis curves of the transitions can be determined. It turns out that the reversal points of the 
control parameter are Markov points which determine the stochastic evolution of the process. It 
has been shown that the branches of the hysteresis loop are converging to fixed limit curves when 
the number of cyclic back-and-forth variations of the control parameter between two consecutive 
reversal points is large enough. This convergence to limit curves gives a clear explanation of the 
accommodation process. The accommodated minor loops show the return-point memory property 
but this property is obviously absent in the case of non-accommodated minor loops which are not 
congruent and generally not closed. In contrast to the traditional Preisach model the reversal 
point susceptibilities are non-zero finite values. The stochastic model can provide a surprisingly 
good approximation of the Raylaigh quadratic law when the external parameter varies between two 
sufficiently small values. The practical benefits of the model can be seen in the numerical analysis of 
the derived equations. On one hand the calculated curves are in good qualitative agreement with the 
experimental observations and on the other hand, the estimation of the joint distribution function 
of the up and down switching fields can be performed by using the measured hysteresis curves. 
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I. INTRODUCTION 

The phenomenon of the hysteresis, understood in a general sense, has been investigated so intensively for many 
decades that any list of references wou^be,-|far from complete by any standards. It is very fortunate that in the 
last few years outstanding monographsEmtrQ have been published in this field, and thus the author does not feel 
obliged to cite the large, amount of old but important references. However, it is considered important to mention 
two papers of G. Kadana'S whose work has played a stimulating role in getting to the idea of reconsideration of thft 
hysteresis theory by the present author. -There is no doubt that the abstract .reformulation of the Preisach modelQ 
given by Krasnoselskii and A. PokrovskiiEl and summarized by I. D. Mayergoyza in his book resulted in an improved 
mathematical clarity in the hysteresis theory, but the stochastic nature of the hysteresis still has not been treated 
with sufficient mathematical rigor. 

The aim of the present paper is to define a stochastic model of hysteresis and to derive exact equations for the 
probability distribution functions describing the state variations in hysteretic processes as a function of increasing as 
well as decreasing control parameters. The vocabulary of magnetic hysteresis will be used for convenience from now 
on, however, the concepts can easily be generalized for any hysteretic phenomenon. 



II. DESCRIPTION OF THE MODEL 



Let us assume that the unit volume of the system consists of many small abstract regions, called "particles" which 
are characterized by four random variables: ju,, A, Xd an d Xu- The absolute value of the particle magnetization is 
denoted by [i. If the particle magnetization is parallel (antiparallel) to the external magnetic field H then the particle 



1 



is in the state (S^) and A = 1(— 1). The random value Xd corresponds to a local field at which the state 
jumps to the state and similarly the =4> transition occurs at the random local field Xu- For simplicity 
the Xu and Xd will be called as U-field and D-field. These two random variables express the obvious fact that each 
particle "feels" not only the external magnetic field H , but also the interaction field due to the adjacent particles and 
the random field originated from the inhomogenities of the surrounding environment. These particles characterized 
by the random variables (i, \,Xd and Xu can be regarded as "independent" abstract elements of the system, and they 
will be called "hysterons". 

Figure [j] illustrates a possible realization of transitions S( — ) of a hysteron. The transition curves form a 

random rectangular hysteresis loop which is almost in all cases asymmetrical in the co-ordinate system of magnetization 
versus external magnetic field since the U- and D-fields are supposed to be random. 

Let us denote by hi, h 2 , ■ ■ ■ , h jy the hysterons in a system, and let be the set of indices and the number 

of hysterons which are in the state <S( + ) at a given external field H. In this case the magnetization of the system is 
given by the stochastic equation 

N 

k=i 

where fik is the absolute value of the magnetization of the hysteron hk, while 

f 1, iffceAA(+), 
A fe = I 

[ —I, if k 3 Af(+l 

Since the random variables (i\, (i%, . . . , (In are mutually independent and have the same probability distribution 
function 

P{(i k < x} = L{x), V k=l,2,...,N, 
it is obvious that the characteristic function of the distribution function 



P{^„(+) < x} = R n (+)(x) 



(1) 



can be written in the form 



* n(+ ,M=/ e iux dR nl+) (x) = [ ( p(-uj)} 



N 



<p(-u) 



,<+) 



where 



+ 00 



e iuJX dL(x) 



+ 00 



e luJX dL{x). 



In order to calculate the characteristic function 

N 



N 



n(+)=0 



n( + )=0 



<p(-u>) 



(2) 



(3) 



(4) 



we need the probability of finding hysterons in the state at the external field H which is the endpoint of a 
well-defined magnetization prehistory. The determination of this probability and the derivation of the equations for 
"up" and "down" magnetizations versus magnetic field will be the task of the next section. 



III. DERIVATION OF THE FUNDAMENTAL EQUATIONS 



A. Some basic relations 



Let us denote by H(x,y\C) the joint distribution function of the random U- and D-fields. From the physical point 
of view it is quite obvious that the U-field cannot be smaller than the D-field, so the stochastic inequality Xu > Xd 
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must be satisfied. It is easy to showu that the joint distribution function of Xu and Xd satisfying the condition 
C = {Xu > Xd} can be written in the form 

£ dx'f* h(x',y') A(x'-y')dy' 
~P{Xu <x, X d< y\Xu > Xd} = H(x, y\C) = °° +QO °° , , (5) 

J-co dx Loo h ( x >V ) d V 

where A(:r) is the unit step function. It is clear that the joint density function of the U- and D-fields can be given by 

h{x,y) 
I-^dx'J^hix'^dy' 



h(x, y\C) = - 2, ^ A(x - y), (6) 



provided that the condition C is valid. 

We need in the sequel two conditional probability distribution functions: 



V{Xu <x\C} = H(x, oo\C) = F u (x\C) = J ; + ~ ; ; J 'f° , ; ,' , , (7) 



f^dx' f^hjx^y^dy' 
dx> r^Hx' ,y>) dy> 

and 



J^dy' J y ^h(x',y')dx' 
f^dy>ti°°h(x',y>)dx>- 



P{Xd < y\C} = ff (oo, y\C) = F d {y\C) = 7 + ~ , , " — . (8) 



Evidently F u (x\C) is the probability that the U-field of a given hysteron is not larger than x, while F^{y\C) is the 
probability that the D-field is not larger than y assuming in both cases that the condition C is fulfilled. By using the 
Dirichlet's theorem for changing the sequence of integration it is obvious that 

+ OC r-x' /- + 00 /- + 00 

dx' h(x',y')dy' = dy' h(x',y')dx'. 



For the sake of further considerations it is necesarry to introduce two transition probabilities denoted by w u (Hi ] H) 
and Wd(H u J, if). Let if; be a fixed value of the external magnetic field and let us suppose that the state of a 
given hysteron is at if. By using elementary theorems, it can be proved that if the external field increases 

monotonically from Hi to H > Hi then 

, fftm i" l dx J-o h (x,y)dy 
w u (Hi t H) = — (9) 

Jh, dx J-ooM^y) d y 

is the probability of the transition iS*- - ' =>• occuring in the interval [Hi | if]. Similarly, let H u be an other fixed 
value and 5 < - + - ) the state of a hysteron at ff u . If the external field decreases now monotonically from H u to if < if u 
then 



w d {H u }H)= J ^ H J » , — (10) 



Ih U d y I y °° h (x,y) dx 
J H ^dy Jy°°h(x,y) dx 

gives the probability of the transition =>■ occuring in the interval [if u J, if] 



B. Stochastic magnetizing process 



Let us introduce the "time parameter" t G [0, +oo] and define a real valued, external field function H(t) which 
consists of monotone increasing and decreasing sections of different length. Denote by ifi, if2, . . . , Hj, . . . , the ex- 
tremum values of the function H(t) belonging to the subsequent time points t\ < t% < ■ ■ ■ < tj < ■ ■ ■. It is clear 
that if H(tj) = Hj is a local maximum then if(t,_i) = if,_i and ff(< J+1 ) = if,+i must be local minimums which 
are not necessarily equal. In the following the sequence {if/} will be called magnetizing path and the elements of 
this sequence are called points of reversal. If the functions H^(t) and H^ 2 '(t) have the same magnetizing path then 
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they are said to be equivalent for any magnetizing process irrespective of the form of the time function between the 
individual extrema. In Fig. ^ two equivalent H(i) functions are seen. The sequence of extrema is the same for both 
curves, but the time distance and the shape of sections between the consecutive extrema are different. 

It is assumed that the magnetizing process which consists of random transitions <^ S^' of hysterons does 
not "feel" the variation speed of H(t) between the consecutive extremum values, i. e., the magnetizing process is 
static. The evolution of the process in each subinterval [Hj, Hj+i], j = 1, 2, . . . , is stochasically determined by the 
extremum Hj and by the actual values of H(t) following Hj, but the process does not depend, on the time derivative 
of H(t). This property is called rate independence in the non-stochastic theory of hysteresisS but it will be applied 
in this stochastic theory too. Denote the maximum reversal fields by odd and the minimum ones by even indices. In 
this case it is clear that 

and naturally any one of the inequalities H 2 k+i > H 2 k-i and H 2k +i < H 2 k-\ can be valid. 

Now, let us define the random function £j + ^(i?2fe-i I H) which gives the number of hysterons in the state at 

the decreasing external field H belonging to the interval [H 2k -\ [H 2k \. Similarly, denote by (H 2 k T H) the number 
of hysterons in the state at the increasing external field H belonging to the interval [H 2k ] H 2 k+\\- 

We suppose that at the starting point of the magnetizing process each hysteron is in the S^ - ' state, that is the 
system is in the state of negative saturation. In the following this fact will be expressed by the stochastic equation 
Quirt = 0. Evidently any other state of the system could as well be chosen for the starting point, this choice, however, 
does not really matter since the influence of the starting state on the evolution of the process — as it will be shown 
— disappears very rapidly. 

In order to describe the magnetizing process we should determine two probabilities. One of them is 

P{^ +) (H 2k+1 i H) = n^ +1 {H)\s}+l t = 0} = p^ +1 [H 2k+1 j H, r$ +1 (H)\0] , (11) 

and the other is 

P{^+\H 2k ^H) = n^(H)\&l t = 0} =p%[H 2k 1Hn%(H)\0}. (12) 
It is important to note that the reversal points (extremum values) 

Hi,H 2 , . . . , H 2 k-i, H 2 k, H 2 k+i, . . . , 
are Markov-points of the stochastic processes 

d +) (Hj t H) and ^ +) {H j+1 J. H), and therefore we can write the following 

equations: 

JV 

Pal+i i H 2k+i I H, n^ +1 (H)\0] = £ p% x [H 2k+1 iH,n^ +l {H)\n^{H 2k+1 )}x 

n l 2 u k } (H 2k + 1 )=0 



xpM[H 2k 1H 2k+1 ,n)%(H 2k+1 )\0}, 



(«)/ 



(13) 



and 



p<$\H 2k \Hnf k \H)n = Yl P 2 ^[H 2 k}H,n 2 u k \H)\n 2 d li(H2k)]x 



I»^_l(H'3t)=0 



x P 2 a U[H2 k -x i H 2k , n^H^O] . (14) 

As it has already been mentioned the hysterons can be regarded as independent of each other particles and, 
therefore, it is an easy task to determine the probability that the number of iSW-hysterons is exactly equal to a non- 
negative integer not larger than N, at an either decreasing or increasing external field H provided that the number 
of iS( + )-hysterons is known at the last reversal point before arriving at H . 

The probability p 2 d J +1 [H 2k +i J, H, n 2 d k+1 (H)\n 2 ^} (H 2k +i)] can be obtained as a result of the following consideration. 
If the number of -hysterons at the reversal point H 2k +\ is equal to n 2 u k (H 2k +\), then — in order to have n^ +1 (H) 
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hysterons in the state at the external field H < i^2fe+i — exactly n^(i?2fe+i) — n 2k+i(H) hysterons of state 
have to transform to the state S^~' in the interval [i?2fc+i [H]. It is obvious that the probability of this event 
can be given by 



piti[H2 k+ iiH >n ^ +1 (H)\n^(H 2k+1 )} = ( n %f 2 , fe +° 

V n 2k+l( H ) 



X 



X [w d (H 2k+1 lH)] ^H^ k+l )-n^(H) {1 Wd{H2k+liH) f£ +i m . (15) 

Similarly, to determine the probability p^} \li2k T H, n^ k (H)\n < ^_ 1 (H2k)) one has to recognize that if the number of 
t 5(-)_hysterons at the reversal point H 2 k is equal to N — n^_ 1 (H2k), then — in order to have n^{H) hysterons 
in the state at the external field H > Hik — exactly n^}(H) — n < ^_ 1 (H2k-i) hysterons of state have to 
transform to the state <S>' + ) in the interval [H2 k T H] . The probability of this event is given by 

P il\H2klH,n^ m nit^k)}=( (X^ti^ 

\ n 2k ( H > - n 2k-l( H 2k) J 

x [w4H 2k lH)] n %W- n %-^ [1 - w u {H2k]H)] N -< )(H) . (16) 
In order to symplify the further calculations let us introduce the generating functions 

N 

T^ +1 (H 2k+1 iH,z) = Ai[H2k + ilH,nW +1 (H)\0] z<trW (17) 

and 

N 

T^{H^H,z)= Y, P { 2 U k\H2k]H^{H)nz^^. (18) 

By using the Eqs. ( |l3| ) and (|l5| ) we get the first fundamental equation in the form 

T^ +1 (H 2k+1 iH,z) = T^[H 2k ^H 2k+1 ,a(H 2k+1 ,H,z)}, (19) 

where 

a(H 2k+1 ,H, z) = w d (H 2 k+ilH) + [1 - w d (H 2k+1 lH)} z. (20) 
The second fundamental equation follows from the relations (jlj) and (|l6|). We have 

T ( ^{H2 k ]H,z) = [c(H 2k ,H,z)] N T ( 2 d k l 1 [H2k-ilH2k,b(H 2 k,H,z)}, (21) 

where 

c{H 2k ,H,z) = l-(l-z)w u {H 2 k]H) 1 (22) 

and 

Now we will derive the characteristic function of the probability that the system magnetization is not larger than 
x at a decreasing external field H which follows the last reversal point H 2 k+i- Introducing the notation 
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and by using the relation ( |l7| ) we obtain from Eq. (0) 

<5> d {H 2k+1 iH,uj) = [<p(-u)] N r^ +1 [H 2k +ilH,^(Lj)]. (25) 

Similarly, if we use Eq. ( |l8| ) then the characteristic function of the probability that the magnetization of the system 
is not larger than x at an increasing external field H after the last reversal point H 2k can be obtained from Eq. (jij) 
in the form 

$„0H2kT-ff,u>) = [<p(-w)] N T { ^[H 2k ]H^{w)]. (26) 

These two characteristic functions describe completely the stochastic behaviour of the magnetizing process in both 
increasing and decreasing external magnetic fields. It is apparent from the above considerations that the stochastic 
model developed by us has been built up without any reference to a particular nature of hysteresis and therefore, its 
generality is at least as high as that of the Krasnoselskii and Pokrovskil3 model. 



C. Calculation of the hysteresis curves 



The expectation value of the magnetic moment /j, k due to the fc-th hysteron can be given by 

d(p(u>) 



E{/i fc } = r 1 



du> 



= M s . 



j=0 



By using this expression we can write the expectation value of the magnetization of the system at the decreasing 
external field H following the reversal point H2k+i m the form 



d$ d (H 2k+1 lH,uj) 



M^ +1 {H 2k+1 [H\Q) 



j=0 



2M S N^ +1 (H 2k+1 iH\0) - N M s , 



where 



Nj£ +1 (H 2k+1 lH\0) 



dT^ +1 (H 2k+1 jH,z) 
dz 



From the fundamental Eq. (O) we obtain 



N^ +1 (H 2k+1 iH\0) ^ N^'(H 2k ^ H 2k+1 \0) [1 - w d (H 2k+l 



(27) 



(28) 



(29) 



The expectation value of the magnetization of the system at the external magnetic field H increasing after the reversal 
field H 2k can be obtained from the equation 



d<£ u (H 2k ^H,Lo) 



dio 



M^{H 2k ]H\Q) 



j=0 



where 



2M S N^{H 2k ]H\0)-N Ms, 



N^(H 2k ^H\0) 



dT^ k \H 2k ]H,z) 
dz 



The recursive relation 

N^{H 2k ]H\Q) = Nw u (H 2k ']H) + N^l 1 (H2k-ilH 2k \0) [1 - w a (H 2k 
follows from the other fundamental Eq. (|2l|). 



(30) 



(31) 



(32) 



G 



By introducing the relative magnetizations 



and 



m^ +1 {H 2k+1 [H\Q) = j^j- M^ +1 (H 2k+1 iH\0), 



(33) 



m 2 u k \H 2k ]H\0) = j^- M^(H 2k ^H\0), 



(34) 



from Eqs. (|27j), (|29j) and (pCf), ( p2[ ) after elementary calculations the following recursive relations are obtained: 

m^ +1 (H 2k+1 iH\0) = 



and 



1 + m^(H 2k ^H 2k+1 \0) [1 - w d {H 2k+1 [H)} - 1, 



m[ u k \H 2k ]H\0) 



(35) 



2w u (H 2k ^H) 



Ad) 
L 2k- 



i^fc-ii^fclO) [l-w^flafcT-ff)]-!- 



(36) 



In order to solve this system of recursive equations we need the formula for the starting branch of the relative 
magnetization. Since the negative saturation has been chosen as the initial state of the system it follows from Eq. (j^) 
that if Hi =>■ Ho = — oo, then 

w u {H^H) => F U {H\C), 
and so the equation for the starting branch will be 



m Q u \-oo]H\0) = 2 F U {H\C) - 1. 



(37) 



This branch can be also called the limiting ascending branch because there is no branch below it. If the positive 
saturation would be the initial state then it is easy to show that the limiting descending branch can be written in the 
form 



m[ d) '(oo I H\0) = 2 F d {H\C)- 1, 



(38) 



where F d (H\C) is defined by the expression (ph, and at the same time it is obvious that the limiting descending branch 
has the property that there is no other branch above it. The two limiting curves form the major hysteresis loop which 
defines an area where all other loops should be located. 

By using the expression ([37j) for the starting branch and Eq. ( |35"1) we can obtain the first descending branch 

mf\H 1 ^H\0) = 2F u {H 1 \C) [1 - w d {H x | H)] - 1, (39) 

which is attached to the limiting ascending branch at the point Hi. This descending branch is called by Mayerg oyzE 
the first-order transition curve. The field Hi where the first-order transition curve starts from, will be called start 
field. Denote by H 2 the next reversal point where the magnetizing field begins again to increase. The corresponding 
ascending branch, i. e., the second-order transition curve is given by the formula: 



>\H 2 ]H\0) = 2w u {H 2 ]H) + [1 + m [ «>(HilH 2 )} [1 - w u {H 2 ] H)] - 1 



(<*)/ 



(40) 



This procedure can be continued and it is seen that there is no need to take into account any special requirement in 
order to describe the field dependence of the average magnetization since the Markov points of the magnetizing field 
determine exaclty the stochastic behaviour of the process. 
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D. Stationarity of hysteresis loops 



Let us investigate now the variation of the magnetization for a special sequences of reversal fields. Let as suppose 
that H2k+i — H u , V k = 0, 1, . . . , while H 2k = Hd, V k = 1,2,..., and H u > Hd, i.e. the magnetizing field is varying 
between two extreme values H u and Hd. The field variation which starts with a decrease of the the external magnetic 
field H from the reversal point H u until it reaches the next reversal point Hd and then turns to increase to the nearest 
H u value, is called magnetizing cycle. The magnetizing cycle results in a hysteresis loop called minor hysteresis loop. 
The first cycle corresponds to the variation of the external field between the reversal points Hi =4> Hi H3 , where 
Hi = H3 = H u and Hi = Hd, while the k-th cycle is done by the variation of the magnetizing field between the 
reversal points H 2k ^i => H 2k => H 2k+ i, where H 2k -i = H 2k+ i = H u and H 2k = H d for k = 1, 2, . . . , . For the 
descending branch of the k-th minor loop one can obtain from Eq. ( ^5|) the following expression: 

m^_ 1 {H u lH\0) = [l + m^_ 2 (H d ^H u \0)} [1 - w d {H u [H)] - 1. (41) 

while for the ascending branch of the k-th minor loop the relation 

m$(H d lH\0) = 2 w u (H d -]H) + [1 + m^^H^lH^O)] [1 - w u (H d ^H)\ - 1 (42) 

can be derived from Eq. (|3(]). By using Eq. (|3^) we have 

m[ d) (H u lH\0) = 2 F V {H U \C) [1 - w d (H u lH)} - 1, (43) 

for the descending branch of the first minor loop and 

m^\H d ^H\0) = 2 F U (H U \C) [1 - w d (H u lH d )} [1 - w u {H d ]H)] + 2 w u {H d ]H) - 1 (44) 

for the ascending branch of the same loop. It is to note that according to the Mayergoyz's terminology the first 
minor loop consists of a first-order descending and a second-order ascending transition curves. Following the reversal 
points of the magnetizing field this procedure can be continued and we can obtain both the descending and ascending 
branches of relative magnetization for any minor loop. 

One can prove a very important limit theorem, namely, there exist two limit curves 

lim m^iHulHlO) = m d (H u lH), (45) 



and 



lim m^{H d ^H\0) = m u {H u {H), (46) 



which are determining a closed minor loop. In other words, the magnetizing process becomes stationary with increasing 
number of cycles. It means the system "forgets" gradually its initial state by repeating the magnetizing cycle. This 
forgetting process can be related to the well-known accommodation process. The original Preisach model results in 
an immediate formation of the minor hysteresis loop after only one cycle of back-and-forth variation of the input 
between any two c onse cutive extremum values. However, this consequence of the Preisach model contradicts to many 
experimental factsl!2l'llJ indicating the accommodation plays an important role in magnetizing processes. In order to 
describe the accommodation process the traditional Preisach model was modified rather artificially in the "moving" 
and the "product" modela!3. In contrast to these the stochastic model contains the phenomenon of accommodation 
inherently and so there is no need for any additional improving the model. 

The formulas for the limit curves defined by Eqs. ( f45| ) and ( f46| ) can be obtained by some elementary calculations 
in following forms: 



and 



where 



m d {H u lH) = 2 Q{H dl H u ) w u (H d -]H u ) [l-w d {H u lH)\ - 1, (47) 
m u {H d ]H) = 1-2 Q(H d ,H u ) w d {H u [Hd) [1 - w u {H d ]H)], (48) 
Q{H dl H u ) = [wd{H u iHd)+w u {Hd]H u )-Wd(H u lHd)w u {H d ]H u )]- 1 . 
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From these Eqs. two important relations can be derived, namely 



m d (H u i H u ) = m u (H d f H u ), and m d (H u [ H d ) = m u {H d f H d ), 

which show that the return-point memory property is fulfiled for the accommodated minor hysteresis loops. It is also 
obvious, that the accommodated minor loops due to the same pair of reversal fields H d and H u > H d are not only 
congruent but identical since the field values H d and H u unambiguously determine the the branches of stationary 
loops. One has to mention that the accommodated branches m d {oo\,H) and m u (—oo'\H) are exactly identical with 
the limiting descending and ascending branches which indicates the consistency of the theory. 

The explicit form of the expressions m d (H u [H) and m u (H d ] H) which describe the descending and the ascending 
branches of the stationary minor loop between two reversal fields H d and H u > H d has a great advantage in numerical 
calculations in comparision with the well-known Everett integral. It is to be noted that the expressions (^?j) and ( |4g| ) 
are suitable to describe not only symmetrical but asymmetrical hysteresis loops too and it is easy to show that 
symmetrical hysteresis loops can be obtained only if the function h(x,y) has a mirror symmetry expressed by Eq. 

h(x >y ) = h(-y,-x). (49) 

In the following the mirror symmetry of h(x,y) will be assumed. 

It is worthwhile to derive the formula for the virgin curve of the magnetization depending on the parameters of the 
density function h(x,y\C). After some simple manipulations we obtain 

™w = ^«) + m2)'-U). m - 1 ' (50) 



where 



si{H) = ^" AX $~°° h{ - XlV) dy (51) 

I-h dx I-oo h ( x >y) d y 



and 



/ m I-* dy C°° h ( X > y ^ dx 
S 2(H) = — (52) 

J-oo d v I y °° K x > y) dx - 



If the function h(x,y) satisfies the symmetry relation (W9) then it is easy to prove that 



and the initial susceptibility defined by Eq. 



lim mn(H) = 



_ dmo(H) _ Jg h(x,0)dx 
X °"-" ** ~ ^ dx J^h^y) dy (53) 

is different from zero in contrary to the classical Preisach model which gives a nonrealistic zero slope of the virgin 
curve at H — 0. 

IV. NUMERICAL CALCULATIONS AND DISCUSSION 

In order to compute the magnetization vs. field curves we have to know the joint density function h(x,y\C) of the 
U- and D-fields. Since these fields are the sum of many small random components it is reasonable to assume that the 
central limit theorem is approximately valid and so the function h(x,y) in h{x,y\C) can be chosen in the form 

h(x, y) 



2tt a 2 yi^Cf 

x exp [(x - H c f -2C r (x- H c ) (y + H c ) + (y + H c f] J , (54) 







where the meaning of the constants H c , a and C r is clear from the elements of the probability theory. Figure |3j 
shows the contour plot of h(x, y\C) defined by (||) for the parameters H c = 0.2, a = 0.6 and C r = 0.5. The contours 
are belonging to the following values of h(x,y\C) = 0.1, 0.2; 0.3(0.05)0.65 and 0.682. The last one is slightly smaller 
than max( I9 ) h(x,y\C) — 0.682923.... The discontinouity along the line y — x — can be clearly seen in the figure. 
In the sequel this formula will be used in all of our numerical calculations provided that the correlation coefficient 
C r is equal to zero, i. e. h{x,y) = f{x)f(—y) where / is the density function of the normal distribution. This Gfse 
corresponds to the product model introduced by G. Biorci and D. PescettiE3 and used consequently by G. Kadan3'Erli3. 

The relative magnetization vs. field curves are shown in Fig. ^. The parameter values used for the calculation 
are H c = 0.4, a = 0.6 and C r = 0. The curves LA and LD correspond to the limiting ascending and descending 
branches, while the curves indexed in the figure by 1, 2, 3, 4 are the first-, second-, third- and fourth- order transition 
curves defined by the reversal points Hi = 1.2, H 2 = —0.8, H 3 = 0.6, H 4 = —0.6. It is worth noting that the all 
information about the past history of the magnetizing process is transfered by the state of the system in the last 
reversal point. For example, the fourth-order transition curve 4 which is plotted in the field interval [—0.6, 1.5], is 
determined by the state in the reversal point H4 = —0.6. 

It is well-known that in the traditional Preisach model the minor loops which describe the cyclic change of the 
magnetization with back-and-forth variation of the magnetizing field between the same two limiting values are con- 
gruent and the formation of a closed minor loop is realized in one cycle, i. e., the accommodation process is absent. In 
contrast to this the stochastic model contains inherently the accommodation process which is clearly demonstrated 
in Fig. 0. For the sake of orientation the limiting ascending branch LA is also plotted in Fig. || where it is seen 
that the descending branch of the first minor loop starts from the point A due to the first reversal field Hi = 0.8 
and after reaching the reversal point H 2 = H& = —0.2 it turns to increase to the point B which corresponds to the 
next reversal field H 3 = H u = 0.8. One can observe that the first minor loop is not closed, the point B where the 
decreasing branch of the second minor loop starts from, occupies a higher position than the point A, and the end 
point C of the increasing branch of the second minor loop is found above the point B but the distance between the 
points C and B is smaller than that between the points B and A. By repeating the magnetizing cycle between the 
reversal fields Hd = —0.2 and H u = 0.8 the difference between the branches of the same type becomes gradually 
negligible, i.e., the branches converge to limit curves which form-finally a closed stationary hysteresis loop denoted by 
LC. The magnetizing curves measured by Carter and RichardsEa on silicon steel (4.3%5'j) are surprisingly similar to 
that plotted in Fig. g. 

In order to demonstrate the speed of the convergence, the non-accommodated relative magnetizations have been 
calculated in the reversal point H u — 0.8 for the subsequent cycles. Figure || shows that the stationary (i. e., the 
limit) value of the magnetization can be very well approached by repeating the cycle 8-9 times in the case of parameter 
values H c = 0.4, a = 0.6 and H u = 0.3. 

The non-accommodated minor loops due to the same pair of reversal fields are evidently not congruent and generally 
are not closed. However, this non-congruency has nothing in common with that introduced and discussed in details 
by KadarBli3. The non-congruency of the non-accommodated minor loops bounded by the same field limits has a 
quite different origin in the stochastic model, namely, the non-equilibrium response of the system for the cyclic back- 
and-forth variation of the external magnetic field between two consecutive reversal points. It is obvious consequence 
of the non-stationarity of minor loops that the return-point memory property is absent in these loops. 

In order to study the properties of non-congruency of this type the first minor loops belonging to different start 
fields are calculated. Denote by AH = H u — Hd the difference between the consecutive reversal fields. For the 
characterization of the non-accommodated first minor loops due to different start fields Hi let us introduce two 
parameters defined by 

W = W(Hi,AH) = max [mi(#i | H) - m 2 (H 1 - AH | H)], 

Hi — i\H ^CH 1 

and 

= 0(Hx, AH) =m 2 (H 1 -AH1H 1 )-m 1 (H 1 lH 1 ). 

The dependence of these parameters on the start field Hi is shown in Fig. ^ for the parameter values: Hd = —0.2, 
AH = 1. The author of the present paper is far not convinced whether the experimental data contradict or support 
the non-congruency of this type because of the lack of careful measurements. 

It seems to be useful to investigate the remanence properties of systems described by the stochastic model. In Fig. 
H the first-order descending curves which start from different points of the ascending limiting branch LA can be seen. 
The curves starting from the points due to the field values H\ = 1.6, H 2 = 1.4, H3 = 1.2, H4 = 1 are plotted to the 
points of remanences Rl, R2, R3, R4 which are obviously different from the stationary (i. e., the accommodated) 
values. The non-accommodated NR and stationary remanences SR versus start field are shown in Fig. 0. As it 
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is seen the non-accommodated remanences can be negative below a critical start field CR since the initial negative 
saturation has a significant effect on the first-order transition curves. The stationary remanence curve SR. calculated 
from the equations p7| ) and ([l8]) is non-negative in all points of the start field interval. 

The influence of the parameter a on the shape of the major hysteresis loop can be seen in Fig. ^[ As it is expected 
the larger is the parameter a the wider is the hysteresis loop, i.e. the larger non-homogeneity in a system (e.g. in a 
magnetic sample) results in a higher " coercive force" . 

It seems to be useful to calculate the accommodated (stationary) hysteresis loops for different pairs of reversal points 
Hd and H u > Hd- The hysteresis loops plotted in Fig. [ll] correspond to the reversal points: Hd = —1.5, H u = 1.5 
(loop ML1), H d = -1, H u = 1 (loop ML2), H d = -0.5, H u = 0.5 (loop ML3). For the calculation we used 
the parameter values: H c — 0.4 and a = 0.6. For the sake of completeness the virgin curve VC calculated by j50| ) 
and the major loop LL bounded by the limiting ascending and descending curves are also shown in the figure. The 
stochastic model clearly shows that all accommodated minor loops corresponding to cyclic inputs between the same 
two consecutive extremum values are not only congruent but simply identical. 

In Fig. |lj three accommodated first-order minor loops denoted by 1, 2, 3 can be seen. The descending branches of 
the loops are started from the field values H — 0.5,0.3,0, and each of the ascending branches returns exactly to the 
same point that the corresponding descending branch left. The returning curves have an apparent slope discontinuity 
with regard to the major loop ALA. 

At this point it is worth to make a remark of somewhat historical nature. As it is well-knpwn Preisach's idea for 
his model was originated from the quadratic Rayleigh relation which can be easily obtaineoo assuming a uniform 
distribution of the U- and D- fields over the "Preisach triangle". It is interesting to note that in the stochastic 
model the calculated hysteresis loops almost perfectly coincide with that calculated by the Rayleigh formula when the 
reversal fields Hd and H u > Hd and so the magnetizing field H £ [Hd, H u ] are sufficiently small. The hysteresis loop 
R defined by reversal points H u = 0.5 and Hd = —0.5 in Fig. [l3| can be very well approximated by the equations 

m d a (0.5 1 H) = C { d) + C{ d) H + C { 2 d) H 2 , 
m u a {Q.h]H) = C£° + C[ u) H + C { 2 U) H\ 

where 

= 0.13035..., 
C[ u) = 0.79301..., 
-C ( 2 u) = -0.54251... 

in the case of parameter values H c — 0.2 and a — 0.6. In Fig. [l3] the squares q correspond to the values calculated 
by the quadratic equations. The excellent aggreement with the curves of the stochastic model indicates that the 
Rayleigh law can be reproduced in a straightforward way in the stochastic model. 

This model differs from the original Preisach model in a very essential point in relation to the reversal point 
susceptibility. Namely, the non-zero initial susceptibility at the turning points is an inherent property of the stochastic 
model, while the traditional Preisach model can produce positive initial slope only if the Preisach function is supposed 
to have a Dirac-delta like singularity n along the boundery of the Preisach triangle, and that is a rather artificial 
requirement introduced by Mayergoyzu. The susceptibility vs. magnetizing field is seen in Fig. |l4| for the parameter 
values: H c = 0.2 and a — 0.6. The shape of the calculated curve can be expected on the basis of physical considerations 
and corresponds to those found experimentally. 

The estimation of the joint density function h{x,y\C) from measured hysteresis curves was beyond the scope of 
our present theoretical consideration. Of course, one may attempt in simple cases to estimate the parameters of a 
plausible density function (e. g., (|54|)) by an appropriate data evaluation procedure. 

V. CONCLUSIONS 

It has been shown that the Preisach model of hysteresis can be replaced by a new model based on exact concepts of 
the probability theory. In this model the phenomenon of hysteresis has been described as a stochastic process defined 
on a set of all possible values of the control parameter the reversal (turning) points of which are Markov points of 
the process. The one dimensional distribution function of the stochastic process has been exactly determined and 
the magnetizations versus up and down magnetic fields have been calculated as expectation values of the stochastic 
process. It has been proven that the magnetizing process becomes stationary with increasing number of magnetizing 
cycles. It means that for the description of the accommodation process there is no need of any artificial auxiliary 
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assumption since the stochastic model contains the phenomenon of accommodation inherently. In general case the 
model is able to describe the symmetric as well as the asymmetric hysteresis. In relatively small magnetizing fields 
the quadratic Rayleigh law can be easily obtained from the equations of the stochastic model. It is important to 
note that the turning point susceptibilities have non-zero finite values in contrary to the traditional Preisach model 
which docs not take consequently into account the random nature of the elementary switching process. Finally, the 
stochastic model shows that all stationary loops corresponding to the same two limiting values of the magnetizing 
field are equivalent but the non-stationary loops arc non-congruent and in general not closed. 

ACKNOWLEDGMENTS 

The author wish to acknowledge the stimulating discussions with Mr. G. Kadar who helped very much in completing 
the present paper. 



3 



1 I. D. Mayergoyz, Mathematical models of hysteresis (Springer- Verlag, Berlin — Heidelberg — New York, 1991) 

2 A. Visintin, Differential models of hysteresis (Springer- Verlag, Berlin — Heidelberg — New Tork, 1994) 
A. Ivanyi, Hysteresis models in electromagnetic computation (Akademiai Kiado, Budapest, 1997) 

4 G. Bertotti, Hysteresis in Magnetism (Academic Press, San Diego, 1998) 

5 G. Kadar, J. Appl. Phys. 61, 4013 (1987). 

6 G. Kadar, Physica Scripta, T25, 161 (1989). 

7 F. Preisach, Z. Phys. 94, 277 (1935). 

8 M. A. Krasnoselskii, A. V. Pokrovskii, Soviet Math. Dokl. 12, 1388, (1971). 

9 L. Pal, Foundation of the probability calculus and statistics (Akademiai Kiado, Budapest, 1995), (in hungarian), Vol.1, p. 122 

10 W. S. Melville, Journ. Inst. Elcctr. Eng. 97(11), 165 (1950). 

11 R. O. Carter, D. L. Richards, Journ. Electr. Eng. 97(11), 199 (1950). 

12 E. Delia Torre, G. Kadar, IEEE Trans. Magn., MAG-23, 2823, (1987). 

13 G. Biorci, D. Pescetti, II Nuovo Cimento, 7, 829, (1958). 

14 R. Becker und W. Doring, Ferromagnetismus (Verlag von Julius Springer, Berlin, 1939), (in german), p. 222 

15 G. Kadar and E. Delia Torre, IEEE Trans. Magn., MAG-23, 2820, (1987). 



FIG. 1. A possible realization of the transition 5* 
FIG. 2. Two equivalent H(t) curves. 

FIG. 3. Contour plot of h(x,y\C) defined by Eq. (6) where h(x,y) is given by Eq. (57) with parameter values: 
H c = 0.2, a = 0.6 and C r = 0.5. 

FIG. 4. Limiting branches LA, LD and magnetization versus field curves starting from the reversal points Hi = 1.2, 
H2 = —0.8, H.i = 0.7 and H4 — —0.6. The magnetization curves are indicated by 1, 2, 3 and 4. 

FIG. 5. Accommodation process of the minor loop in consecutive magnetizing cycles between the reversal points H d = —0.2 
and H u = 0.8. The loop LC is the accommodated minor loop. 

FIG. 6. Convergence of the relative magnetization in the reversal point H u — 0.8 with increasing number of cycles to the 
limit (i. e., the stationary) value. 
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FIG. 7. Width W of the first-order minor loops and the difference O between the values of the descending and ascending 
branches in the reversal point H u — 0.8 versus start field. 

FIG. 8. The limiting ascending branch LA and four first-order descending curves ending in non-accommodated remanences 
denoted by Rl, R2, R3 and R4. 

FIG. 9. The non-accommodated NR and the stationary SR relative remanences versus start field due to different points of 
the ascending limiting branch. 

FIG. 10. Influence of the parameter a on the shape of the major hysteresis loop in the case of H c = 0.4. 

FIG. 11. The virgin curve VC, the major loop LL and three accommodated hysteresis loops ML1, ML2, ML3 calculated 
for different pairs of reversal points in the case of parameter values H c = 0.4 and a = 0.6. 



FIG. 12. Three accommodated first-order minor loops denoted by 1, 2, 3 and the shifted ascending branch ALA. 



FIG. 13. The hysteresis loop between "small" reversal points and the quadratic Rayleigh curves denoted by □. 



FIG. 14. The irreversible susceptibility versus magnetizing field H. 
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